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Abstract 


Return periods and the probability of the occurrence of precipitation are two important indicators 
in studies related to flooding and precipitation, and a correct understanding of the indicators can 
play an important role in climate risk studies. This study aimed to investigate the performance of 
the Multi-Source Weighted-Ensemble Precipitation (MSWEP) and to use it in determining the 
return periods and the probability of annual precipitation in Iran. For this purpose, the generalized 
extreme value (GEV) method was used. Four statistics of RMSE, MBE, PBIAS, and R? in 
Köppen-Geiger climate classification in Iran were used to evaluate the MS WEP data. The results 
showed that this data has the best performance in arid and semi-arid regions of the country with 
a PBIAS of 0.40 and 0.32. In contrast, this data showed the least performance in two rainy (Cfa) 
and mountainous areas (Dsb), so that the maximum PBIAS between the climate zones of the 
country is seen between -2.45 and -0.09% in the above-mentioned zones, respectively. The results 
showed that the maximum probable precipitation occurs in periods of 1 to 15 years in the northern 
coasts and its maximum intensity occurs in the southern coasts of Iran. The probability of 
occurrence of precipitation in Iran fluctuates between 2.03 to 32.02%. Decreasing latitude from 
north to south and decreasing height from west to east are associated with a decrease in the 
probability of occurrence of precipitation in Iran. Maximum uncertainty is seen in the 
mountainous regions of Iran with a high probability of occurrence of precipitation. 
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جغرافیا و مخاطرات محیطی, سال دهم شمارة چهلم. زستان ۱۶۰۰ صص ۲۰۹-۲۲۷ 


مقاله پژوهشی 


Å + ty 
MA GALEN 


بررسی دوره باز گشت 3 احتمال رخداد بارش در ایران بر اساس پایگاه داده بارش همادی -وزنی جند 
منبعی (MSWEP)‏ 


I‏ زرین - استادیار اقلیم‌شناسی, گروه جغرافیاء دانشگاه فردوسی مشهد 
عباسعلی داداشی رودباری- پژوهشگر پسادکتری اقلیم‌شناسی. گروه جغرافی؛ دانشگاه فردوسی مشهد 


تاریخ دریافت: ۱2۰۰/۲/۳۱ تاریخ بازنگری: ۱۶۰۱/۵/۱۰ تاریخ تصویب: ۱۶۰۰/۷/۱۹ 


حکیده 


دوره بازگشت و احتمال رخداد بارش دو نمایه مهم در مطالعات مربوط به سیل و بارش‌های فرین 
محسوب می‌شوند و درک صحبح از این دو نمایه می تواند نقش مهمی در مطالعات ریسک اقلیمی داشته 
باشد. این پژوهش با هدف بررسی کارایی پایگاه داده بارش همادیسوزنی چند منبعی (MSWEP")‏ و 
کاربست آن در تعیین دوره‌های بازگشت و بررسی احتمال رخداد بارش سالانه در ایران انجام شده است. 
برای این منظور از روش توزیع فرین تعمیم يافته GEV")‏ استفاده شد. برای درستی سنجی داده‌های 
پایگاه 1۷15۱۷7۳ از چهار سنجه آماری PBIAS MBE RMSE‏ و R?‏ در پهنه‌های اقلیمی کوپن- 
گایگر در ایران استفاده گردید. نتایج نشان داد که این پایگاه داده در مناطق خشک و نیمه‌حشک کشور با 
درصد اریبی ۰/٤١‏ و ۰/۳۲ بهترین کارایی را دارد. در مقابل» این پایگاه داده در دو پهنه پرباران (Cfa)‏ و 
کوهستانی (Db)‏ کمترین کارابی را خود OLE‏ داد؛ به‌طوری‌که بیشینه درصد اریبی در بین پهنه‌های اقلیمی 
کشور به ترتیب با ۲/۶۵-و tA‏ درصد در پهنه‌های مذ کور دیده شد. نتایج نشان داد که بیشینه مقدار 
بارش محتمل در دوره‌های ۱ تا ۱۵ ساله در سواحل شمالی و بیشینه شدت آن در سواحل جنوبی ایران 
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2 Multi-Source Weighted-Ensemble Precipitation (MSWEP) 
3 Generalized Extreme Value (GEV) 
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11۰ جغرافیا و مخاطرات محبطی شمارة چهلم 


افا ی او ال خرن ی e Ea‏ ت ی کا ان راا ی وا 
در ایران در ارتباط هستند. بیشینه عدم قطعیت در احتمال رخداد بارش روزانه در مناطق کوهستانی ایران 
مشاهده شد. 


کلیدواژه‌ها: ایران» پایگاه داده بارش MSWEP‏ توزیع فرین تعمیم (GEV) asl‏ دوره بازگشت بارش. 


۱- مقدمه 


مخاطرات طبیعی تهدید مهمی برای امنیت عمومی. زیرساخت‌هاء منابع طبیعی و توسعه اقتصادی در سراسر جهان 
هستند (راگنو و همکاران؛ ۲۰۱۹). در سال‌های اخیر» رخدادهای فرین اقلیمی " در بسیاری از نقاط جهان به میزان 
قابل توجهی افزایش lash‏ (الکساندر و همکاران " ۲۰۰7؛ وال و همکاران" ۲۰۱۵؛ مایره و همکاران“ NG‏ سال و 
Coll See‏ ۲۰۲۱). به flo‏ اهمیت فرین‌های اقلیمی در مدیریت ریسک و طراحی مهندسی (جنگ و همکاران " 
۶ علاقه زیادی به شناخحت چگونگی تغییر در این رخدادها به‌ویژه دوره بازگشت بارش و احتمال رخداد OF‏ 
وجود دارد. بررسی دوره‌های بازگشت بارش و احتمال رخداد آن برای کشوری همچون ایران با اقلیم غالب خشک و 
نیمه‌حشک می‌تواند برای مدیریت منابع آب و پیش‌بینی احتمال وقوع سیلاب بسیار مفید باشد. 

از روش‌های آماری ایستا و نا ایستا برای بررسی ارتباط رخدادهای فرین بارشی (به‌عنوان‌مثال بارش‌های سنگین و 
خیلی سنگین و سیل) و احتمال وقوع آن‌ها استفاده می‌شود (استدینگر و همکاران/ ۱۹۹۳). در مطالعات اقلیمی؛ 
تکنیک‌های تحلیل فراوانی کلاسیک و ایستا نتایج درستی را از واقعیت‌های اقلیمی ارائه نمی‌دهند. چراکه احتمال 
رخداد بارش و دوره بازگشت از رویکرد ناایستایی پیروی می‌کنند (خالق و همکاران“ ۲۰۰7). بااین‌حال, فقط تعداد 
کمی از مطالعات بارش را با مدل‌های ناایستا در ایران بررسی کرده‌اند. ازاین‌رو لزوم یک مطالعه جامع در این حوزه 
eet oe‏ 

افزایش رخدادهای فرین و پیامدهای گزارش شده از آن‌ها در سراسر جهان انگیزه استفاده از روش‌های نایستا را 
برای اطمینان از به دست آوردن مشخصههای در حال تغییر اقلیم افزایش داده است (میلی و همکاران "4 ۲۰۰۸). 


بااین‌حال» نظرات و دیدگاه‌های گوناگونی در مورد مدل‌های tlb‏ برای تحلیل فرین‌ها وجود دارد. بسیاری از 
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روش‌های موجود برای تحلیل فرین‌ها که تحت فرض نایستایی می‌باشند. دارای محدودیت‌های زیادی ازجمله فقدان 
یک چارچوب کلی برای ترکیب متغیرهای درگیر در برآورد پارامترهاو عدم قطعیت در برآورد احتمال رخداد 
فرین‌های بارشی و دوره‌های بازگشت وجود دارد (گائو و همکاران ۲۰۱۹). بررسی دوره‌های بازگشت و احتمال 
رخداد بارش برای طیف گسترده‌ای از کاربردهای پژوهشی. مهندسی و عملیاتی ضروری است. از اطلاعات دوره‌های 
بازگشت بارش در بخش پژوهشی می‌توان برای ارزیابی خروجی مدل‌های اقلیمی» در بخش مهندسی برای 
برنامه‌ریزی و در زمینه عملیاتی نیز برای طراحی و مدیریت زیرساخت‌ها استفاده کرد. همچنین برای بهبود سامانه‌های 
پیش‌بینی سیل و حشکسالی نیز می‌توان از این رویکرد مطالعاتی استفاده نمود (بوسیلوویج و همکاران » ۲۰۰۸؛ هیرپا 
و همکاران ؛ ۲۰۱7؛ بان و همکاران " ۲۰۱۹). توانایی برآورد قابل‌اطمینان از مقدار و فراوانی مورد انتظار رخدادهای 
فرین برای ome‏ عرصه‌های مدیریتی ریسک» اساسی است. از طرف دیگر برآورد دوره باز گشت و احتمال رخداد 
بارش از موضوعات مهم در مطالعات اقلیمی است. بر این اساس ضمن استفاده از روش‌های کارآمد برای بررسی 
رخداد فرین‌های بارشی, استفاده از داده‌های معتبر نیز یک ضرورت جدی تلقی می‌شود. پایگاه‌های داده بارش دارای 
اریبی زیادی در مناطق کوهستانی هستند و ازآنجایی که بخش‌های بزرگی از ایران در شسمال‌غرب و غرب را مناطق 
مرتفع کوهستانی تشسکیل می‌دهند. استفاده از داده‌ای که بتواند کمترین خطا را ارائه دهد از اهمیت شایان توجهی 
برخوردار اسست. داده‌های همادی از Glee‏ کمتری در برآورد بارش نسبت به داده‌های منفرد برخوردارند (بیک و 
همکاران, ۲۰۱۹) بر این اساس با تکیه بر نتایج داده‌های همادی می‌توان نتایج مطمئن‌تری را برای برآورد احتمال 
رخداد بارش و دوره بازگشست OT‏ به‌دست آورد. دوره‌های بازگشست بارش و احتمال رخداد آن یکی از موضوعات 
موردتوجه در علوم جوی و بخصوص اقلیم‌شناسی و هیدرولوژی است. همان‌طور که پیش‌تر بحث شد دو رویکرد 
ایستایی و اایستایی در بررسی فرین‌های اقلیمی موردتوجه است. نظر به اهمیت رویکرد ناایستایی و استفاده از این 
روش در پژوهش حاضر در ادامه پژوهش‌هایی که از رویکرد gl wll‏ استفاده کرده‌اند موردتوجه و بررسی قرار 
می‌گیرند. در یکی از این پژوهش‌هاء گرویزمن و همکاران * (۲۰۰۵) با بررسی روند بارش فرین در اقلیم‌های مختلف 
کره زمین بر اساس روش GEV‏ نشان دادند که مناطقی با مقدار بارش زیاد در دوره‌های بازگشت. دارای رخدادهای 
بارش سنگین و بسیار سنگین هستند. در پژوهشی دیگر دگانتو* (۲۰۰۹) دوره‌های بازگشت بارش در ایالات‌متحده 
آمریکا را با روش GEV‏ برای دوره آماری ۱۹۵۰-۲۰۰۷ موردبررسی قرار داد و به این نتیجه رسید که امکان دارد 
طوفان‌های بزرگ در هر *۶ سال یکبار اتفاق بیافتند. همچنین ایشان نشان داد که بارش‌هایی با دوره باز کشت ۲ ساله 
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به ازای هر دهه ۲ درصند tol fl‏ یافت‌اند. گائو و همکاران ۱ (۲۰۱۳) a‏ مدل‌سازی ab ull‏ بارش‌های فرین در چین 
با استفاده از 1۳۱ ایستگاه طی دوره ۱۹۵۱-۲۰۱۳ پرداختند و به این نتیجه رسیدند که توزیع‌های GEV‏ نایستا کارایی 
بهتری از معادل ایستایی of‏ دارند. برآورد پارامترهای GEV‏ و تحلیل ایستایی و ناایستایی از بارش‌های فرین در آفریقا 
پژوهشی است که دپائولا و همکاران " (۲۰۱۸) بدان پرداختند. نتایج آن‌ها نشان داد که روند بارش از رویکرد نایستایی 
پیروی می‌کند و بررسی اثرات تغییر اقلیم با رویکرد ناایستایی نتایج بهتری را نسبت به روش ایستایی ارائه می‌دهد. وو 
و میشرا" (۲۰۱۹) با بررسی فراوانی رخدادهای فرین بارش در آمریکا ضمن آشکارسازی افزایش فراوانی رخدادهای 
بارش فرین نان دادند که بیش‌تر بارش‌های فرین از حالت ناایستایی پیروی می‌کنند. ونر و همکاران * (۲۰۲۰) به 
برزسی ذوره‌های باز گشت فرین‌هاق دما و بازش بر اساس پرونداد مدل‌های 01۷11۳6 پرداعتند و به این جه زسیدند 
که با افزایش طول دوره بازگشست فراوانی فرین‌های اقلیمی کمتر می‌شسود. بارش و پایگاه‌های داده بارش موردتوجه 
پژوهش‌های بسیاری در ايران بوده است. به‌طوری‌که پژوهش‌های بسیاری به‌درستی سنجی پایگاه‌های داده بارش 
و 0۳00 (میری و همکاران؛ ۱۳۹۵ ERA-interim‏ (رضینئی و (YAT «agi‏ و ERA-Interim‏ 
PERSIANN-CDR‏ و CHIRPS‏ (گرجی زاده و همکاران ۱۳۹۸) پرداختند. بااین‌حال تا کنون پژوهشی که داده‌های 
همادی را در سطح کشور درستی سنجی و استفاده کرده باشد. انجام نشده است. همان‌طور که تاکنون پژوهشی از 
قابلیت‌های داده‌های وزنی چند گانه برای برآورد بارش در کشور استفاده نکرده است» کمتر پژوهشی نیز روش‌های 
ناایستایی را برای بررسی بارش و رخدادهای فرین آن موردتوجه قرار داده است. از معدود پژوهش‌هایی که با استفاده 
از روش فرین تعمیم GEV) th‏ به بررسی بارش‌های فرین در ایران پرداخته‌اند می‌توان به پژوهش عساکره (۱۳۹۱) 
و مظفری و همکاران (NAD)‏ اشساره کرد. عساکره (۱۳۹۱) تغییر توزیع فراوانی بارش‌های فرین شسهر زنجان را با 
استفاده از توزیع فرین تعمیم يافته GEV)‏ موردبررسی قرار داد. نتایج وی نشان داد که چند نوع اایستایی در 
بارش‌های زنجان رخ داده اسست. این ناایستایی‌ها به صورت چرخه‌ها؛ جهش و نیز تغیبر توزیع فراوانی قابل‌رژیت 
است. مظفری و همکاران OTAN‏ نیز با واکاوی روابط فضایی بارش‌های فرین در غرب Ol pl‏ نشان دادند که الگوی 
فضایی بارش‌های فرین در غرب ایران یک الگوی خحوشه‌ای است که معنی‌داری آن در سطح ٩‏ درصد تأیید 
می گردد. 

جمع‌بندی پیشینه پژوهش نشان می‌دهد که به‌رغم عدم قطعیت پایین‌تر روش‌های ناایستایی و درک مؤثرتری که 
اا وی ی داشاو کا ی ی وی کی کیش Gs‏ ایران ریخات ا فرار_کرفته نیش 


این پژوهش دوره از ets‏ و احتمال رخداد آن در ol yl‏ را بر اساس روش‌های ناایستایی با کاریست داده‌های پایگاه 
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پژوهش با سه هدف اصلی ۱- درستی سنجی داده‌های پایگاه بارش همادی-وزنی MSWEP‏ در مناطق مختلف 
اقلیمی ایران؛ ۲-بررسی دوره بازگشت و ۳- بررسی احتمال رخداد بارش با استفاده از روش توزیع فرین تعمیم asl‏ 


توجهی برخوردار است. 


۲- مواد و روش 
aisge‏ 


ایران کشوری وسیع در جنوب غرب آسیا است و به جهت اقلیمی کشوری خشک و نیمه‌خشک محسوب 
می‌شود. ارتفاع ایران از کمتر از VA‏ متر در سواحل جنوبی دریای زر تا OWS‏ متر در قله دماوند در تغییر اسست 
(شکل ۱). رشته‌کوه‌های البرز در شمال و زاگرس در غرب نقش بسزایی در تقسیم کشور به پهنه‌های مختلف اقلیمی 
دارند. این مناطق کوهستانی همانند سدی در مقابل رسیدن رطوبت به قسمت مرکزی کشور عمل می NS‏ که درنتیجه 
بارش در این مناطق نسبت به سایر مناطق کشور بسیار کم‌تر است. ایران همچنین میزبان یکی از گرم‌ترین بیابان‌های 
جهان-دشت لوت- است. بر اساس روش طبقه‌بندی اقلیمی کوپن بیش از ۸۰ درصد از مساحت کشور در مناطق 
Sis‏ و نیمهخشک واقع شده‌اند (رضیئی» AFA‏ 
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شکل ۱- طبقه‌بندی اقلیمی ایران بر اساس روش کوپن -گایگر (تفکیک افقی ۱ کیلومتر) (بیک و همکاران 
۸ و مدل رقومی ارتفاعی (SRTM)‏ با تفکیک افقی ۲۰ متر به همراه پراکنش ۱۲۰ ایستگاه همدید طی دوره 


آماری ۱۹۹۰-۲۰۱۷ 


1٤‏ جغرافیا و مخاطرات محبطی شمارة چهلم 


(MSWEP) پایگاه داده بارش همادی-وزنی چند منبعی‎ -Y-Y 


در بین متغیرهای هواشناسی» بارش به دلیل ناهمگنی زمانی - مکانی یکی از پیچیده‌ترین متغیرها است (پرین و 
گوبیت» ۲۰۱۷). داده‌های مشاهداتی به‌عنوان منبع اصلی جهت بررسی بارش‌های فرین هستند؛ بااین‌حال تعداد 
ایستگاه‌های هواشناسی در ایران که داده بلندمدت بارش در اختیار داشته باشند بسیار محدود است و در کویرهای 
bel‏ و مناطق کوهستانی ایران عملاً ایستگاه‌هایی با داده بارش بلندمدت وجود ندارد. طی دهه‌های گذشته مجموعه 
داده‌های بازتحلیل بسیاری همانند 0۳0 MERRA-2 CRU" ERAS"‏ و داده‌های ماهواره‌ای همانند 112۷1۷ 
GPM‏ و PERSIANN-CDRY‏ برای بارش ایجاد شده‌اند (سان و همکاران. ۲۰۱۸). این مجموعه داده‌ها از نظر منبع 
داده (ایستگاه رادار زمینی ماهواره بازتحلیل شده و یا همادی» تفکیک افقی (از ۰/۰۵ تا ۲/۵ درجه قوسی) گام 
زمانی (۳۰ دقيقه تا ماهانه) و از نظر هدف از بازتولید مجموعه داده‌ها (دقت مکانی. همگنی زمانی» طول دوره یا 
ترکیبات (Of‏ متفاوت هستند. هر داده شبکه‌ای بارش حتی داده‌های ماهواره‌ای دارای عدم قطعیت و بیش‌برآوردی /کم 
برآوردی هستند. بهترین راهکار برای کاهش عدم قطعیت داده‌های شبکه‌ای استفاده از الگوریتم همادی‌سوزنی است؛ 
چراکه تلفیق داده‌های مختلف باعث کاهش اریبی داده‌ها می‌گردد (ببشاپ و آبراموویتز“ ۲۰۱۳). لذا به این دلیل در 
این پژوهش از داده‌های بارش پایگاه همادی-وزنی MSWEP‏ استفاده گردید. پایگاه داده بارش همادی-وزنی چند 
منبعی یک مجموعه داده جهانی بارش است که دارای گام‌های زمانی ۳ ساعته. روزانه و ماهانه است. در این پژوهش 
از 1۷151۷۳-۷220 استفاده گردید که دارای دوره آماری ۱۹۷۹ تا ۲۰۱۷ است Se)‏ و همکاران» ۲۰۱۹). این 
مجموعه داده نسبت به ساير داده‌های شبکه‌ای بارش ازاین‌جهت منحصربه‌فرد است که از یک رویکرد ترکیبی -وزنی 
مبتنی بر داده‌های ایستگاهی, ماهواره‌ای و بازتحلیل بهره می‌برد تا برآورد قابل اطمینان‌تری از بارش را در کره زمین 
ارائه دهد. از زمان انتشار نسخه ۱ با تفکیک افقی ۰/۲۵ درجه قوسی) در ماه می ۲۰۱۳ L MSWEP‏ موفقیت در 
مقیاس جهانی برای اهداف مختلفی مورداستفاده قرار گرفته است. نسخه دوم پایگاه داده MSWEP‏ تلفیقی از داده‌های 
مبتنی بر ایستگاه (GPCC  WorldClim, GHCN-D,GSOD)‏ ماهواره CMORPH, GridSat, GSMaP, TRMM)‏ 
و (TMPA3B42RT‏ و بازتحلیل ERA-Interim)‏ و JRA-55‏ است. پایگاه داده بارش MSWEP-V220‏ از چندین 


جنبه نسبت به سار داده‌های موجود منحصربه‌فرد است: ۱) تفکیک افقی (۰/۱ درحه قوسی) و گام زمانی الا (۳ 


1 Global Precipitation Climatology Centre (GPCC) 

2 ECMWF Reanalysis ۷5 (ERAS) 

3 Climatic Research Unit (CRU) 

4 Modern-Era Retrospective analysis for Research and Applications, version 2 (MERRA-2) 

5 Tropical Rainfall Measuring Mission (TRMM) 

6 Global Precipitation Measurement (GPM) 

7 Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks - Climate 
Data Record (PERSIANN-CDR) 

8 Bishop & Abramowitz 

9 Beck et al 
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ساعته)؛ ۲) تلفیق وزنی طیف گسترده‌ای از مجموعه داده‌های ایستگاه. ماهواره و بازتحلیل‌هاء برای برآورد بهتر بارش 
در هر مکانء ۳) تصحیح اریبی برای باران‌ریزه (Drizzle)‏ ۶) تصحیح اریبی سیستماتیک بارش با ایستگاه‌های زمینی 
به دلیل کمبود ایستگاه در ۱۳۷۹۲ حوضه آبریز در سراسر جهان؛ ۵) تصحیح اریبی بارش با استفاده از مشاهدات 
روزانه از ۷۱۷۶۷ ایستگاه در سراسر کره زمین (بیک و همکاران ۲۰۱۹). 

۳-۲- ایستگاه‌های همدید ایران 

برای درستی سنجی برونداد بارش پایگاه ۷15۷۷۳۳-۷220 از داده‌های ایستگاه‌های هواشناسی همدید استفاده 
شد. برای گزینش ایستگاه‌های همدید سه شرط موردتوجه قرار گرفت: ۱-مشارکت حداکثری ایستگاه‌ها در Wl pl‏ 
-Y‏ حداکثر داده مفقود از کل سری زمانی ایستگاه بیش از ۱۰ درصد نباشد (زولینا؛ ۲۰۰۳)؛ و ۳- همگنی داده‌ها با 
آزمون الکساندرسون (SNHT)‏ (الکساندرسون" OAT‏ بررسی و همگنی آن‌ها در سطح ۰/۰۵ مورد تأیید قرار گیرد. 
این آزمون ازآن‌جهت برای همگنی داده‌ها مورداس‌تفاده قرار گرفت که آزمونی با بازدهی بالاست که می‌تواند به دو 
شکل مطلق و نسبی مورداستفاده قرار گیرد (رحیم‌زاده و نساجی‌زواره» ۱۳۹۳). در این پژوهش برای بررسی همگنی 
سری‌های سالانه بارش از نوع مطلق این آزمون استفاده گردید. بر اساس سه شرط یاد شده ۱۲۰ ایستگاه برای دوره 


آماری ۱۹۹۰-۷ برای درستی سنجی بارش گزینش گردید که در شکل (۱) موقعیت آن‌ها نشان داده شده است. 
¥-2— سنحه‌های آماری مورداستفاده جهت درستی سنجی 


در این پژوهش از چهار سنجه آماری PBIAS® MBE! RMSE"‏ و 22 استفاده شد که جزئیات آن‌ها در جدول 


seal مه‎ OD 
حدول ۱- سنحه‌های آماری مورداستفاده در این پژوهش‎ 
نحوه محاسبه آماره‎ RTE 
1 * 
RMSE = ور‎ ۷2-۱0۶ Ki) = ZK) (RMSE) مجذور میانگین مربعات خطا‎ 
n 
n — 
MBE = 2-00 — P) (MBE) میانگین اریبی خطا‎ 
n 
n obs _ sim 
Pbias = =f a (PBIAS) درصد مقدار اریی‎ 
i=1( i 
BE CRO EEE 
È(P; — P)2(0; — 0)? A 


1 Zolina 

2 Alexandersson 

3 Root Mean Square Error (RMSE) 
4 Mean Bias Error (MBE) 

5 Percentage of Bias (PBIAS) 

6 Coefficient of Determination (R2) 
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در روابط MBE‏ و P; «R‏ مقدار برآورد ان مه Oj‏ مقدار اندازه‌گیری» و 0 میانگین مقادیر اندازه‌گیری شده 
می‌باشند. همچنین در RMSE‏ (:6) 77 مقدار برآورد شده متغیر در L(X;) Xi‏ مقدار مشاهده‌ای متغیر در an Xi‏ 
تعداد کل داده‌های مشاهده‌ای است. در رابطه Yoim (PBIAS)‏ مقدار برآورد شده (پایگاه (MSWEP‏ و Yops‏ مقدار 
مشاهداتی (ایستگاه‌های همدید منتخب) است. 

۵-۲- دوره‌های بازگشت بارش در ایران با روش poles‏ فرین تعمیم (GEV') ash,‏ 

بررسی فراوانی بارش‌های فرین و دوره بازگشت آن‌ها از اهمیت بالایی در مدیریت ریسک سیلاب و خشکسالی 
برخوردار اسست SL)‏ و همکاران ‏ ۵۹ که می‌توان آن را با استفاده از توزی‌های احتمال با نظریه مقدار فرین 
برآورد کرد (کولز " ۲۰۰۱). نظریه مقدار فرین EVT)‏ به‌عنوان یکی از شاخه‌های مهم علم آمان به‌طور گسترده در 
علوم هواشناسی و هیدرولوژی مورداستفاده قرار گرفته است (سوکیسیان و تسالیس* ۲۰۱۵). در بیش‌تر موارد. EVT‏ 
اعمال شده در هواشناسی و هیدرولوژی بر این فرض استوار است که این رخدادها ایستا" هستند (کاتز و همکاران " 
۲ بااین‌حال» فرض ایستایی به دلیل بروز تغییر اقلیم و وردایی آن در سال‌های اخیر به‌تدریج به چالش کشیده 
شده aal‏ (سالاس و (Yeg ral Swal‏ ناایستایی ممکن است بر شدت و فراوانی رخدادهای فرین تأثیر بگذارد 
AY)‏ ۲۰۱۳؛ مونیر و گائو OT‏ ۵ بنابراین» مدل‌های توزیع احتمال ناایستایی '' بايد موردبررسی فراز گیرند. دو 
روش فرین تعمیم (GEV) 8b‏ و توزیع پارتو تعمیم یافته (GPD)‏ به‌طور گسترده‌ای برای تحلیل رخدادهای فرین 
مورداس‌تفاده قرار می گیرند (آکرو و همکاران ؛ ۰۱۰۰ رولفوا و همکاران ۲ ۰۳۰۹+ يوزگاتلیگیل و ترکش ۸ T1۸‏ 
SE‏ و همکاران* ۲۰۱۹؛ خاویر و همکاران" ؛ ۲۰۱۲۰؛ کیم و همکاران"/ ۲۰۲۰). در این مطالعه برای محاسبه دوره 


بازگشت مجموع بارش سالانه در öll‏ از روش فرین تعمیم (GEV) abl‏ استفاده شد. 


1 Generalized Extreme Value (GEV) 
2 Fang et al 

3 Coles et al 

4 Extreme value theory (EVT) 
5 Soukissian and Tsalis 

6 stationary 

7 Katz et al 

8 Salas and Obeysekera 

9 Mika 

10 Monier and Gao 

11 Nonstationary 

12 Acero et al 

13 Rulfova et al 

14 Yozgatlıgil & Türkeş 

15 Fang et al 

16 Xavier et al 

17 Kim et al 
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دوره‌های نار گس به‌عنوان aly‏ کمی نوزیع GEV‏ محاسبه pel tye ae‏ تابع Je‏ تجمعی ) (CDF'‏ (آقاکوچک و 
همکاران, ۲۰۱۱؛ خارین و همکاران ؛ ۳ توزیع GEV‏ همانند معادله (۱) است. 


)\( 0 عد € 


GEV (x) = exp |- ca + gaa 


سطح بازگشت Xp‏ مربوط به دوره‌های بازگشت L)‏ بازه‌های بازگشت) Vy‏ = 7 را می‌توان برآورد SoS‏ در آن 


۳- یافته‌ها 

۱-۳- درستی سنجی پایگاه داده بارش MSWEP-V2‏ در ایران 

برای درستی سنجی برونداد بارش پایگاه MSWEP‏ در مقیاس سالانه از چهار سنجه آماری MBE RMSE R?‏ 
و PBIAS‏ استفاده شد. همچنین برای بررسی بهتر کارایی این پایگاه در مناطق مختلف اقلیمی. متوسط هر یک از 
آماره‌ها در پهنه‌های مختلف اقلیمی کوین-گایگر ارائه شد. روش پهنه‌بندی اقلیمی کوپن-گایگر به دلیل مقبولیت بالا 
در سطح جهانی نسبت به ساير روش‌های موجود طبقه‌بندی اقلیمی از یکسو و همچنین در نظر گرفتن متغیرهای 
به این دلیل است که دوره بازگشت و احتمال رخداد بارش در مقیاس سالانه محاسبه می گردد که بر این اساس بهتر 
می‌توان نسبت به نتایج به‌دسست‌آمده در مقیاس سالانه بحث کرد. مقدار ضریب تبیین R’)‏ برای برونداد بارش 
MSWEP‏ در ایران بین ۰/۸۳۲ تا ۰/۹۹ در تغییر است که نشان‌دهنده کارایی قابل‌قبول این پایگاه داده بارش در کشور 
است. بالاترین ضریب تبیین در مناطق حشک و نیمه‌حشک ايران به‌دست آمده است که با توجه به مقدار کم بارش 
سالانه در این مناطق جندان دور از انتظار نبوده استت: همان‌طور که در شکل o> (Y)‏ می‌شود بيشینه مقدار RMSE‏ 
RMSE‏ نیز با مقدار کمتر از ۱ میلی‌متر در جنوب‌شرق و مناطق داخلی ایران دیده می‌شود. با توجه به نتایج بررسی 
این دو سنجه می‌توان نتیجه گرفت که این پایگاه داده در مناطق خحشک و نیمه‌عشک کارایی بسیار بالاتری دارد. 


1 Cumulative Density Function (CDF) 
2 Kharin et al 
3 Kharin and Zwi 
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مقادیر MBE‏ نیز همانند RMSE‏ است و نشانگر آن است که از جنوب به شمال Ol yl‏ مقدار اریبی افزایش می‌یابد. 
در مقیاس سالانه اریبی بارش در ایران بین WEE‏ تا ۱۵/۷۱ میلی‌متر در نوسان اسست. در زاگرس مرتفع و جنوب 
غربی دریای خزر مقدار اریبی منفی شده است (کم برآوردی). بر این اساس می‌توان اذعان داشت که MSWEP‏ در 
مناطق پربارش و مرتفع ایران دارای کم‌برآوردی اسست. از ۱۲۰ ایستگاه موردبررسی برای درستی سنجی پایگاه 
در ۵7/۶۱ درصد از ایستگاه‌ها اریبی منفی (کم‌برآوردی) و در ۹ درصد از ایستگاه‌ها اریبی مثبت 
(بیش‌برآوردی) به‌دست آمده است. به‌طورکلی بررسی اریبی در پهنه‌های اقلیمی نشان می‌دهد که بارش در پهنه‌های 
Dsa Csa BWk BSh BSK Cfa‏ و Dsa‏ دارای اریبی منفی است (جدول (Y‏ که نشان دهنده کم برآوردی این 
پایگاه بارش در پهنه‌های مذکور است. 

همان‌طور که در شکل (۲) و جدول (Y)‏ مشخص است بر اساس سنجه‌های MBE RMSE‏ بارش در مناطق 
شمالی Cha)‏ و (Csa‏ و زاگرس مرتفع Dsb)‏ دارای اریبی و خطای زیادی است که مبین کارایی نه‌چندان مناسب این 
پایگاه داده در مناطق شمالی ایران پا به عبارتی در پهنه‌های اقلیمی Csa Cfa‏ و Dsa‏ است؛ اما مقایسه نتایج این دو 
سنجه آماری با ۳131۸5 نتیجه متفاوتی را نشان می‌دهد. علت این تفاوت در نحوه محاسبه دو سنجه RMSE‏ و MBE‏ 
است. ازآنجاکه RMSE polis‏ و MBE‏ تابعی از دامنه تغییرات بارش در هر ایستگاه اسست. مقدار زیاد RMSE‏ و 
MBE‏ در ایستگاه‌های ناحیه خزری (پهنه‌های اقلیمی Cha‏ و (Csa‏ واقع در شمال ايران (ایستگاه‌های بندر cd pl‏ 
بابلسر گرگان» رامسر و رشت) به جهت مقدار زیاد بارش در این ایستگاه‌ها است. لذا به کاربستن تنها دو سنجه 
MBE‏ و RMSE‏ نتایج صحیحی را ارائه نمی‌دهد. برای رفع این مشکل از سنجه PBIAS‏ استفاده شد. 

نتایج مربوط به سنجه PBIAS‏ نشان داد که برغم مقدار بالای MBE‏ برای پهنه‌های شمالی (مقدار MBE‏ به‌دست 
آمده در پهنه اقلیمی Cfa‏ ۵۸/۳ - میلی‌متر است) این مقدار نسبت به کل بارش این پهنه‌ها کمتر از ٤‏ درصد (درصد 
اریبی منفی) است که نشانگر دقت مناسب داده‌های ۷15۷۷۳۳-۷220 در پهنه‌های خزری Ol yl‏ است. بر عکس مقدار 
MBE‏ که تابعی از دامنه تغییرات بارش است و با افزایش بارش مقدار آن به شکل غیرواقعی افزايش می‌یابد» لذا 
5 سنجه کاراتری برای برآورد درصد اریبی بارش در ایران است. از نظر پراکنش فضایی درصد اریبی (PBIAS)‏ 
بارش در چهار منطقه از Ol pl‏ منفی است: ۱- سواحل جنوبی دریای خزر؛ ۲- رشته‌کوه‌های جنوبی البرن ۳- 
ارتفاعات زاگرس و ۶- جنوب شرقی ایران. بیشینه درصد اریبی منفی بارش در ارتفاعات زاگرس قابل‌مشاهده است. 
در مقابل در مناطق شمال‌غربی. شمال‌شرقی و مناطق بزرگی از ایران داحلی درصد اریبی MSWEP‏ مثبت است. 
به‌طورکلی کوهرنگ با درصد اریبی “VN‏ درصد بيشینه درصد اریبی منفی و بم با ۷/۹۳ درصد بیشینه درصد اریبی 
مثبت را نشان داده‌اند. همان‌طور که روشن است اگر تنها به دو سنجه RMSE‏ و MBE‏ بسنده می‌شد به نظر 
می‌رسید که بم نسبت به بندر انزلی اریبی کمتری دارد؛ اما همان‌طور که مشخص شد مقدار درصد اریبی از کل بارش 


سالانه در بندر انزلی تنها ۳/۷- درصد است که به‌رغم کم‌برآوردی در بارش دارای دقت بهتری نسبت به بم است. 


شال دهم بررسی دوره ناز گشت و احتمال رخداد بارش در asss‏ ۳۹ 


جدول ۲- متوسط سنجه‌های آماری در پهنه‌های اقلیمی ایران بر اساس روش کوپن-گایگر 
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شکل ۲- توزیع فضایی سنجه‌های MBE RMSE R2‏ و ۳11۸5 بر اساس مقایسه بارش سالانه پایگاه داده 


MSWEP-V2‏ و بارش ایستگاهی در ایران 


۳۳۰ جغرافیا و مخاطرات محبطی شمارة چهلم 


-Y-Y‏ دوره بازگشت مجموع بارش سالانه در ایران 

دوره بازگشت بارش سللانه و احتمال رخداد بارش‌های برابر و بیش از ۱ میلی‌متر در ایران با استفاده از روش 
GEV‏ محاسبه و نتایج آن در شکل (۳) ارائه شد. بر اساس روش GEV‏ بارش سالانه ایران دارای دوره‌های بازگشت؛ 
۵ ۳ ۵ ۰ و ۱۵ ساله است. بیشینه شدت بارش در دوره‌های بازگشت اشاره شده در سواحل جنوبی کشور زاگرس 
مرتفع و سواحل خزری دیده می‌شود. 

از نظر پارامترهای GEV‏ افزایش مقدار بارش با افزایش گام زمانی دوره‌های بازگشت در مناطق جنوبی ایران 
(مناطق ساحلی خلیج‌فارس و دریای عمان) همزمان با افزایش پارامتر مکانی توزیع 017۷ است (شکل C‏ بر این 
اساس می‌توان اذعان داشت که روند افزایشی در دوره‌های بازگشت در سطح آماری ۵ درصد )0.05 = 4) معنی‌دار 
است. 

بیشینه مقدار بارش در دوره‌های بازگشت ۱ ۳ ۵ ۱۰ و ۱۵ ساله افزایش زیادی در سواحل خلیج‌فارس و دریای 
عمان نشان می‌دهد (شکل ۳). در این مناطق مقدار بارش از دوره بازگشت ۱ ساله به ۲ ساله با Kal‏ متوسط ۶ درصد 
به ازای هرسال افزایش می‌یابد. مطابق با شکل (۳) افزایش مقدار بیشینه بارش با افزایش دوره بازگشت در سواحل 
جنوبی ایران بارزتر است اما بايد توجه داشت که درست در جنوب‌غربی ایران و مناطق شمالی استان خوزستان مقدار 
بارش در دوره‌های بازگشت موردبررسی به ازای هر دهه حدود ۱ درصد کاهش می‌پابد. 

همان‌طور که در شکل (۳) مشخص است و پیش‌تر نیز گفته شد با افزايش دوره بازگشت. بیشینه مقدار بارش نیز 
افزایش می‌یابد. این افزايش در سواحل جنوبی ایران کمی بیشتر از مناطق شمالی کشور است؛ اما به‌طورکلی در بین 
مناطق و دوره‌های مختلف آهنگ افزایش بارش با طولانی‌تر شدن دوره‌های بازگشت همراه است. 

بیشینه مقدار بارش در دوره بازگشت ۱ ساله در Ol pl‏ بین ۷/۹۵ میلی‌متر در جنوب شرق تا ۸۲/۰۳ میلی‌متر در 
زاگرس مرتفع در نوسان است. با افزایش دوره ES GL‏ مقدار بارش محتمل در یک روز نیز افزايش می‌یابد. 
به‌طوری که بیشینه بارش در دوره بازگشت ۳ ساله ۱۱۳/4۹ میلی‌مت ۵ ساله ۱۲۸/۲۸ میلی‌مت ۱۰ ساله ۱۵۳/۳۰ 
میلی‌متر و ۱۵ ساله ۱۱۰/۸ میلی‌متر برآورد شد. احتمال بارش در دوره‌های بازگشت به‌دست آمده در سواحل جنوبی 
ایران از چابهار تا اهواز, زاگرس مرتفع در غرب ایران (استان‌های لرستان» کرمانشاه و ایلام) و سواحل خزری نشان از 
خطر بالای سیل در این مناطق از ایران دارد. در سال‌های pH!‏ سیل‌های ویرانگری در مازندران (۱۳۹۷/۱۲/۲۳ بر 
اساس گزارش سازمان هواشناسی کشور طی دو روز ۳۰۰ میلی‌متر بارش داشته است که باعث کشته شدن ۱۳ نفر 
شد و به تعطیلی فرودگاه» راه‌آهن و تخلیه چند روستا انجامید)» گیلان (۸۱۳۹۹/۰۵/۱۵ ۳۷۸ میلی‌متر بارش در ۲۶ 
ساعت و خسارت به ۶۰ تا 10 هزار هکتار از زمین‌های کشاورزی) پل‌دختر (۰۱۳۹۸/۰۱/۰۵ مطابق با گزارش سازمان 
هواشناسی. ۱۰۸ میلی بارش در یک روز, ۱۵ کشته و خسارت به ۳۳۰۰ خانه) و سیستان و بلوچستان (۱۳۹۹/۰۱/۰۶ 


مطابق با گزارش سازمان هواشناسی ٩‏ میلی‌متر بارش در یک روز در نیکشهر) اتفاق افتاده انیت مقایسه سیلاب‌های 


شال دهم بررسی دوره باز کشت و احتمال رخداد بارش در NTA asss‏ 


رخ داده با مقدار بارش در دوره بازگشت‌های مختلف که با روش GEV‏ محاسبه (Ad‏ نشان از دقت بالای این روش 
در محاسبه دوره‌های باز گشت بارش در ایران دارد (زشکل ۳. بر این اساس مناطقی را که با رنگ قرمز در نقشه نشان 


داده clo‏ می‌توان پهنه‌هایی h> L‏ سیل‌خیزی بالا عنوان 55 
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در ایران برای مجموع بارش سالانه‎ (GEV) شکل ۳- دوره بازگشت بارش با استفاده از روش فرین تعمیم یافته‎ 


مبتنی بر برونداد پایگاه MSWEP‏ 


۳-۳- احتمال رخداد بارش و عدم قطعیت OT‏ در Ol gl‏ 

احتمال رخداد بارش یک میلی‌متر و بیش‌تر معباری برای شناسایی شدت بارش‌های روزانه است. بالا بودن 
شاخص احتمال رخداد بارش و پایین بودن عدم قطعیت. نشان از احتمال بالای بارش در یک روز را دارد. تأثیر 
عرض جغرافیایی و ارتفاع بر احتمال رخداد بارش‌های روزانه در اران مشهود است. نتایج این نشان داد که با کاهش 
عرض جغرافیایی از شمال به جنوب مقدار این نسبت کاهش می‌یابد. استثنائاتی نیز دیده می‌شسود که مربوط به 
ازتفاعات داشعلی کشسور است که باعت افزایش استمال رخداد بارش می‌شود. در همین راستا خوالفقاری و همکاران 


۳۳۲ جغرافیا و مخاطرات محبطی شمارة چهلم 


(۱۳۸۸) نشان دادند که عرض جغرافیایی و ارتفاع رابطه معنی‌داری با حداکثر بارش‌های روزانه و بارش‌های سالانه 
در ایران دارد. 

عدم قطعیت احتمال رخداد بارش روزانه در ایران در تمامی مناطق کوهستانی بالای ۰/۵ است. نتایج به‌دست آمده 
پیش‌تر در بخش درستی سنجی نیز مورد تأیید قرار گرفته است. به‌طوری‌که پهنه اقلیمی 250 با 8/94 درصد اریبی, 
بالاترین متوسط اریبی را در بین بهنه‌های اقلیمی Ol pl‏ داشسته است. بر این اساس می‌توان اذغان داست که افزایش 
ارتفاع در ایران سبب کاهش قطعیت در برآورد احتمال رخداد بارش روزانه می‌شود. 

seer ده ان‎ AUC راهن آیران‎ E که تیال تفارش ها یک سای‎ Cie, alt 
است که پراکندگی درصد احتمال رخداد بارش با توزیع فضایی بلندمدت بارش در ایران کاملاً همخوانی دارد.‎ 
به‌طوری که بیشینه درصد احتمال رخداد بارش با بیش از ۳۰ درصد در سواحل جنوبی دریای خزر و کمینه آن با‎ 
حداکثر ۲ درصد در جنوب شرق ایران مشاهده می‌گردد.‎ 

به جهت جغرافیایی چهار منطقه را می‌توان به‌عنوان مناطقی با احتمال بالای رخداد بارش روزانه در ایران در نظر 
گرفت. این چهار منطقه عبارتند از: ۱- سواحل جنوبی دریای خزر, ۲- مناطق مرزی غربی ایران از آذربایجان غربی تا 
کردستان. ۳- زاگرس مرتفع در کوهرنگ بختیاری و -٤‏ مناطق مرتفع شمال‌شرقی ایران در مرز با کشور ترکمنستان. 
در مقابل جنوب شرقی و سواحل OLS‏ به همراه بخش‌های بزرگی از ایران مرکزی پایین‌ترین احتمال رخداد بارش 
روزانه را دارند. در این مناطق احتمال رخداد بارش یک میلی‌متر و بیش‌تر در یک روز کمتر از ۵ درصد است که 
علت آن دوری از منابع رطوبتی و قرارگیری در سایه بارش زاگرس است. 

همان‌طور که در شکل (4) دیده می‌شود این پایگاه داده بر روی بخش‌های مرتفع زاگرس درصد احتمال پایینی را 
برای بارش روزانه برآورد کرده است که [Jo‏ اصلی آن به‌طوری‌که در بخش درستی سنجی داده‌ها نیز اشاره شده 
است کمبود اطلاعات مربوط به ایستگاه‌های زمینی در دوره موردمطالعه در این ناحبه است. 

بر اساس نتایج به‌دست آمده احتمال رخداد بارش یک میلی‌متر و بیش‌تر از OT‏ در ۷۶/۷۱ درصد از کشور کم‌تر از 
Vs‏ دو زور ell‏ له galls‏ که اال تاه بارشی ور Coals SN olf he Ol‏ در بلاق از lle‏ 
۵ درجه شمالی قرار دارند. همچنین مناطق مرتفع زاگرس نیز شامل این مناطق می‌باشند (شکل .)٤‏ 

یافته‌های این تحقیق OLE‏ می‌دهد که عدم قطعیت احتمال رخداد بارش در مناطق مرتفع کوهستانی ایران بالاست. 
به‌طوری که در این مناطق عدم قطعیت تا ۰/۱ نیز می‌رسد. در مقابل مناطق خحشک داخلی و شرقی ایران پایین‌ترین 
عدم قطعیت را با مقدار ۰/۰۷ نشان دادند (شکل .)٤‏ 


شال دهم بررسی دوره از کستت 3 احتمال رخداد بارش در .... ۳۳۳ 
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شکل ۶- احتمال رخداد بارش یک میلی‌متر و بیش تر و عدم قطعیت OT‏ در ایران برای مجموع بارش سالانه مبتنی 
بر برونداد پایگاه MSWEP‏ 


۶- نتیجه گیری 

تعیین مقادیر و فراوانی رخدادهای بارش فرین نقش کلیدی در برنامه‌ریزی مدیریت ریسک و مهندسی دارد. با 
توسعه الگوریتم‌های برآورد بارش و افزایش تعداد پایگاه‌های در دسترس در مقیاس جهانی با تفکیک افقی بالاء امکان 
برآورد صحیح بارش در مناطق مختلف بخصوص مناطق بیابانی و مرتفع کوهستانی فراهم شده است. این تحقیق با 
هدف درستی سنجی و بررسی احتمال دوره‌های بازگشت بارش با استفاده از داده‌های پایگاه ۷15۷۷۳۳ انجام شده 
است. از توزیع فرین تعمیم یافته (GEV)‏ برای ارائه درکی مناسب از الگوهای فضایی دوره بازگشت و احتمال 
رخداد بارش ۱ میلیمتر و بیش‌تر در ایران استفاده گردید. 

پایگاه داده MSWEP‏ که کارایی آن در ایران موردبررسی قرار گرفته اسست. مبین کم‌برآوردی بارش در بیش‌تر 
پهنه‌های اقلیمی کشسور اسست. به‌طوری‌که از ۸ پهنه اقلیمی اصلی ایران که اطلاعات ایستگاه‌های هواشناسی آن در 
دسترس بوده است اریبی بارش به‌طور متوسط در ۷ پهنه منفی (کم‌برآوردی) و فقط در پهنه BSK‏ دارای متوسط اریبی 
مثبت (بیش‌برآوردی) است. بررسی درصد اریبی بارش نسبت به مجموع بارش سالانه در ایران نشان از کارایی قابل 
قبول این داده در سطح کشور دارد. به‌طوری‌که درصد اریبی بارش در ۵ پهنه از ۸ پهنه موردبررسی کمتر از ۱ درصد 


مجموع بارش سالانه است. بیش‌ترین اختلاف بین داده‌های مشاهداتی (ایستگاه‌های همدید) با داده‌های پایگاه 


۳۳ جغرافیا و مخاطرات محبطی شمارة چهلم 


MSWEP-V220‏ در مناطق مرتفع زاگرس دیده می‌شود. این اختلاف ممکن است به این دلیل باشد که این پایگاه از 
داده‌های ایستگاه مبنا WorldClim, GHCN-D,GSOD‏ و GPCC‏ استفاده می‌کند و به دلیل اینکه ایستگاه‌های زاگرس 
مرتفع تازه تأسیس هستند و این چهار پایگاه داده -ایستگاه مبنا- مشارکت يافته در پایگاه ۷15۱۷3۳-۷220 از این 
داده‌ها به دلیل کوتاه‌مدت بودن استفاده نکرده‌اند. 

با توجه به Sul‏ درمجموع پایگاه داده بارش ۷15۱۷۳۳-۷220 کارایی قابل قبولی را در پهنه کشسور ایران از خود 
نشان داد؛ از این داده جهت محاسبه دوره بازگشت بارش و احتمال رخداد بارش‌های یک میلی‌متر و بیش‌تر استفاده 
گردید. نتایج نشان داد که بارش در اران دارای دوره‌های بازگشست ۱ call ۵ calle ۳ calle‏ ۱۰ ساله و ۱۵ ساله است. 
بیشینه بارش محتمل در یک روز به‌ترتیب ۸۲/۰۳ ۱۱۳/6۹ ۸۱۲۸/۲۸ ۱۵۳/۳۰ و ۱۷۰/۸۲ میلی‌متر برای دوره‌های 
باز گشست مذکور می‌باشد که این بیشینه‌ها در مناطق شمال ایران. سواحل خلیج‌فارس و مناطق جنوب‌غربی ایران 
مشاهده می گردند. این مناطق مبین پهنه‌های بارشی سنگین و خیلی سنگین در ایران نیز هستند. در همین راستاه 
گرویزمن و همکاران" (۲۰۰۵) نشان دادند که مناطقی با مقدار بارش زیاد در دوره‌های بازگشست. بارش سنگین و 
بسیار سنگین را تجربه می‌کنند. عدم قطعیت احتمال رخداد بارش روزانه در ایران از حداقل کمتر از ۱ درصد در 
مناطق خحشک داخلی تا بیش از ۰ درصد در مناطق مرتفع کوهستانی Oly!‏ در تغییر اسست. احتمال رخداد بارش 
روزانه با افزایش عرض جغرافیایی و ارتفاع در ایران افزایش می‌یابد. تأثیر عرض جغرافیایی و ارتفاع در پیشینه رخداد 
بارش پیش‌تر توسط ذوالفقاری و همکاران (۱۳۸۸) در ارتباط با نسبت بیشینه بارش‌های روزانه به بارش‌های سالانه 
در ایران مورد تأید قرار گرفته اسست. همچنین عدم قطعیت بالا در مناطق مرتفع کوهستانی به دلیل نبود ایستگاه‌های 
مشاهداتی در پایگاه‌های ایستگاه مبنای مشارکت atl‏ در پایگاه MSWEP-V220‏ چندان دور از انتظار نبوده است. 

تغییرات مکانی احتمال رخداد بارش روزانه نه‌تنها می‌تواند بر مدیریت منابع آب و خطر طغیان سیل و خحشکسالی 
تأثیر بگذار بلکه در تغییر توزیع بارش در طول سری زمانی نیز دارای اهمیت است. اگرچه وردایی زیادی در بیشینه 
مقدار بارش در دوره‌های بازگشت ۱ تا ۱۵ ساله ایران مشاهده می‌شود اما ماندگاری کانون‌هایی با شدت بارش بالا 
قابل‌توجه است که نشان دهنده افزايش خطر بروز سیلاب به دلیل احتمال بیشتر رخدادهای بارش سنگین و خیلی 
سنگین می‌باشد. همچنین ارتباط معنی‌داری بین مناطقی با شدت بارش بالا در دوره‌های بازگشت مختلف و پهنه‌هایی 
با سیل‌خیزی بالا در ایران (واقفی و همکاران. ۲۰۱۹) دیده می‌شود. بر این اساس مناطقی که در این پژوهش Ol ped‏ 
کانون‌های بیشینه بارش در دوره‌های بازگشت شناسایی شدند را می‌توان به‌عنوان مناطقی با خطر سیل‌خیزی بالا در 
نظر گرفت. 


1 Groisman et al 
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